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Abstract. We use low-degree acoustic modes obtained by the BiSON to estimate the main- 
sequence age f Q of the Sun. The calibration is accomplished by linearizing the deviations from 
a standard solar model the seismic frequencies of which are close to those of the Sun. Formally, we 
obtain the preliminary value f = 4.68 ± 0.02 Gy, coupled with an initial heavy-element abundance 
Z = 0.0169 ± 0.0005. The quoted standard errors, which are not independent, are upper bounds 
implied under the assumption that the standard errors in the observed frequencies are independent. 
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1. INTRODUCTION 

Seismological calibration of stellar models against observed frequencies of low 



degree mod e s was first di s cussed more than two decades ago (|Christensen-Dalsgaard 



1984 1 19881 : lUlrichl Il986l : iGoughl Il987h . a nd can be regarded as a means of deter- 



mining the main-sequence age of the Sun (iGuenthed 119891: iGoueh & Novotnvll l990: 
Guent her & Dem arauelll997l: IWeiss & Schlattllll998l : bzie mbows ki et al.lll999T lGough 
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200 ll : iBonanno. Schlattl & Paternq 120021) . The procedure is to match certain appro- 



priate seismic signatures of theoretical frequencies determined on a grid of stellar 
models with corresponding signatures obtained from the observations. The signatures 
are chosen to reflect principally the properties of the energy-generating core, where 
nuclear transmutation leaves behind an augmenting concentration of helium, lowering 
the sound speed relative to the environs and thereby providing a diagnostic of age. 
But the signatures are also susceptible to other properties of the stellar interior, which 
must be eliminated before a robust outcome can be achieved. For example, although 
the so-called small frequency separation is sensitive predominantly to the evolving 
stratification of the core, its dependence on the zero-age chemical abundances plays a 
significant contaminating role. Therefore it behoves us to seek an additional diagnostic 
to attempt to measure abundance separately. For given relative abundances of the heavy 
elements, the total absolute heavy-element abundance Z and the 4 He abundance Y 
are related by the requirement that the model has the observed luminosity and radius, 
principally the former. Therefore we need to aim at detecting only one of them. Here 
we use a signature indicative of the abrupt variation of the first adiabatic exponent J\ 
induced by the ionization of helium. 




FIGURE 1. Top left: The symbols (with error bars obtained under the assumption that the raw frequency 
errors are independent) represent second differences, A2V, of low-degree solar frequencies from BiSON. 
Top right: The symbols are second differences A2 V of adiabatic pulsation eigenfrequencies of solar Model 
S. The solid curve in both panels is the diagnostic (f2j) - ©, whose eleven parameters have been adjusted 
to fit the data optimally. Bottom: The symbols denote contributions 8v to the frequencies produced by the 
acoustic glitches of the Sun (left panel) and Model S (right panel). 



2. THE SEISMIC DIAGNOSTIC AND CALIBRATION METHOD 

Any abrupt variation in the stratification of a star (relative to the scale of the inverse 
radial wavenumber of a seismic mode of oscillation), which here we call an acoustic 
glitch, induces an oscillatory component in the spacing of the cyclic eigenfrequencies 
V n> i of seismic modes, of order n and degree /. Our interest is principally in the glitch 
caused by the depression in the first adiabatic exponent ji = (d\np/dlnp) s (where p, 
p and s are pressure, density and specific entropy) caused by helium ionization. The 
deviation 

ov:=v-v s (1) 

(from now on we omit the subscripts n, I) of the eigenfrequency from the corresponding 
frequency v s of a similar smoothly stratified star is the indicator of Y that we use in 
conjunction with the indicators of core structure to determine the main-sequence age. 

Approximate expressions for the frequency contrib utions 8v arising from acoustic 
glitches in solar-type stars were recently presented by iHoudek & Goughl (|2007|) . Here 
we improve them by adopting the appropriate Airy functions Ai(— x) that are used as 
comparison functions in the JWKB approximations to the oscillation eigenf unctions. 
The complete expression for 8v is then given by 



<5v = 8 r v + 8 c v , 



(2) 





A 


C -Sfi/Yi 


Sun (BiSON) 


0.2764 


1.785 0.04325 


Model S 


0.2780 


1.818 0.04511 



TABLE 1. Asymptotic fitting coefficients A, C 
(see equation|6]l and —8y\/y\ = An/V27TVoAn. 



where 



8 r v = -v^AnAjj 1 [v + i(m + l)v ] 
L Jo 

+ / r K- n 1 e-( T - T ") 2 / 2A n|^| 1 / 2 |Ai(-x)| 2 dT 



X 



(3) 



arises from the variation in fi induced by helium ionization, and 

8 c v ~ A c vlv- 2 {l + l/l6n 2 xlv 2 y l/2 

x |cos[2vA c +tan" 1 (4^T v)]-(16^ 2 f c 2 v 2 +l) 1/2 } 



(4) 



results from the acoustic glitch at the base of the convection zone. Here, m = 3.5 is a 
constant, being a representative polytropic index in the expression for the approximate 
effective phase \y appearing in the argument of the Airy function, and /3, r\ and /! are 
constants of order unity which account for the relation bet ween the acoustic glitche s 
caused by the first and second stages of ionization of helium ( Houdek & Goughll2007h : 
T is acoustic depth beneath the seismic surface of the star, and T = 1 /2vo is the total 
acoustic radius of the star; X\\ and An are respectively the centre and the width of 
the Hell acoustic glitch. The argument of the Airy function is x = sun(i//)|3i///2| 2 / 3 , 
where i//(t) = kcot - (m+ l)cos _1 [(m+ l)/<of] if f > T t , and i//(t) = \k\cot - (m + 
l)ln[(m+ \)/cof + \k\] if t < T t , in which f = T + co^En, with co = 2nv, and T t is the 
location of the upper turning point; k(t) = [1 — (m+ l) 2 /(0 2 f 2 ] 1//2 , and Kj = K"(r/Tn) and 
Kji = k-(th)- Also ty/ c = K c coi c - (m+ l)cos _1 [(m+ l)/f c co] + 7t/4, where K c = k(t c ) 
and f c = T c + CO l £ c . The seven coefficients An, An, Tn, £n, A c , T c , £ c are found by fitting 
the second difference 

3 

A 2 v = V n _ 1;/ - 2v nJ + Vn+ij ~ A 2 (5 r v + 8 c v) + £ a k v~ k (5) 

yt=0 

to the corresponding observations under the assumption that the errors in the frequency 
data are independent (see top panels of Figure 1). The last term in equation© approx- 
imates smooth contributions arising, in part, from wave refraction in the stellar core, 
from hydrogen ionization and from the superadiabaticity of the upper boundary layer of 
the convection zone, introducing four more fittin g coefficients (k =0,...,3). 

The outcome of the fitting to the BiSON data (Bas u et al.ll2007h and to the adia bati- 
cally computed eigenfrequencies of solar Model S (Christensen-D alsgaard et aDl 19961) is 
displayed in Figure 1 : the upper panels display the second differences, together with the 
fitted formula ©, the lower panels display the corresponding contributions 8v to the fre- 
quencies of oscillation from the acoustic glitches. To the resulting smooth ( glitch free) 
frequencies V s , derived from equation (1), of both the solar observations (|Basu et al 




FIGURE 2. Functional forms fx of the integrands <px 
in X = / ^xdr, where fx (r) = 0x ('") / 1 0x (0) | and where 
X = A, C, or F, plotted over the inner half of the inter- 
val (0,R) of r. The parameters A, C and F are sensitive 
particularly to the structure of the core, being progres- 
sively more centrally concentrated. For the calibrations 
here we use only A and C (for F is more difficult to de- 
termine with confidence), whose values determined from 
the fit are listed in Table 1, together with the implied max- 
imum depression — 8ji/Y\ = An/\/2jrvoAii in y x caused 
by He II ionization. 



2007) and the eigenfrequencies of the reference solar model (Model S) was fitted the 



asymptotic expression 

, A , AL 2 -B , CL 4 -DL 2 +E d FL 6 -GL 4 +HL 2 -I fi 
v„,z~(n+iZ + e)v v£ - 3 v£ v o> ( 6 ) 

where L 2 = /(/ + 1), from which we obtain the coefficients Vo, £, A,B,C,D,E,F,G,H 
and / , each of which is an integral of a function of the equilibrium stratification, some of 
which are displayed in Figure 2. The differences between the actual smoothed frequen- 
cies v s and the asymptotic expression © are plotted in Figure 3. 
We have carried out age calibrations using combinations of the parameters 

Za = (A,C,-8YL/YL), a = 1,2,3. (7) 

Presuming, as is normal, that Model S is parametrically close to the Sun, we consider 
the solar value £ ® to be approximated by a two-term Taylor expansion of £ a about the 
value for Models S: 

where Af and AZ are the deviations of age f and initial heavy-element abundance Z 
from ModelS, and e^ a are the formal errors in the calibration parameters. A (paramet- 
rically local) maximum-likelihood fit (again, assuming that the errors in the observed 
frequencies are independent) then leads to the following set of linear equations: 

H uj C al H fik®0k = H ajC a pA Q p , (9) 

in which 0^ = (A? ,AZ) + £©£ = 0o£ + £©/b k= 1,2 is the solution vector subject to 
(correlated) errors £©£, A^ = — ^ + = Ao/3 + C a p is the covariance matrix 

of the errors £^ a , mdH aj = [(d£ a /dt) z , (d£ a /dZ) t ], j = 1,2. 
A similar set of equations is obtained for the formal errors £@£ : 

H aj C ap H pk £ ek = H aj C a p£^, (10) 



from which the error covariance matrix C©^ = £®k £ ®q can be computed with a Monte 
Carlo simulation. 
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FIGURE 3. Differences between the smoothed frequencies v s of the Sun (left) and Model S (right), and 
the fitted asymptotic expression ©. Modes of like degree / are connected by solid lines. 



TABLE 2. Partial derivatives H a j obtained from two sets of calibrated evolutionary models for the 
Sun. Values with respect to age t® are in units of Gy . 



(dA/dt & ) z 


(dA/dZ) te 


(dC/dt @ ) z 


(dC/dZ) t@ 


[d(-8n/Yi)/dt @ ]z [d(-8n/Yi)/dZ} te 


-0.0469 


-0.584 


0.677 


36.8 


-0.00656 0.442 



The partial derivatives H a j are obtained from two sets of five calibrated evo- 
lutionary models for the S un, computed with the evolutionary programme by 



Christensen-Dalsgaardl (|1982t) . and adopting the Livermoore equation of state and 
the OPAL92 opacities. One set of models has a constant value for the heavy-element 
abundance Z = 0.02 but varying age; the other has constant age but varying Z. The 
values of the partial derivatives H a j are listed in Table 2. 

3. RESULTS 

Age calibrations using the different combinations of the parameters t, a are summarized 
in Table 3; error contours associated with the first entry are plotted in Figure 4. In all 
cases the age found is greater than currently accepted values. The values of Z should not 
be regarded strictly as statements about the initial he avy-element abundan ce, but rather 



as measures of the opacity in the radiative interior. lAsplund et al.l (|2004) have argued 
that the photospheric abundances of C, N and O had previously been overestimated, 
suggesting that the actual total heavy-element abundance is rather lower than previously 
believed. However, that cannot imply that the opacity in the solar interior is necessarily 
comparably lower because it has been implicitly calibrated here (by accepting the tenets 
of solar-evolution theory, and the OPAL opacity calculations upon which the models 
are based), and indeed the opacity has already been determi ned seismolo gically from 



a broader spectrum of modes than has been adopted here (|Goug h 2004). The matter 
raised by Asplund etal. therefore challenges either the opacity calculations, the nuclear 
reaction rates, or the basic physics of stellar evolution, not helioseismology, as some 
spectators have surmised. As we know already from seismological structure inversions, 
the solar models are not accurate by helioseismological standards. Therefore the prop- 



TABLE 3. Age calibrations with different combinations of |, 





to (Gy) 


r l/2 

'-on 


Z 


r l/2 
*-022 


-(-Cei 2 ) 1/2 


A,C,-8ft/Yi 


4.679 


0.017 


0.0169 


0.0005 


-0.0023 


A,C 


4.658 


0.023 


0.0177 


0.0007 


-0.0037 


A, -571/71 


4.673 


0.017 


0.0165 


0.0007 


-0.0019 


C,-8yi/ri 


4.700 


0.028 


0.0169 


0.0005 


-0.0029 



0.018 r 



0.017 



0.016 




4.62 4.64 4.66 4.68 4.70 4.72 4.74 
to (Gy) 



FIGURE 4. Error ellipses for the calibration using all 
three parameters E, a : solutions (f©,Z) satisfying the fre- 
quency data within 1, 2 and 3 standard errors in those 
data reside in the inner, intermediate and outer ellipses, 
respectively. 



erties inferred from these calibrations could be more contaminated by systematic error 
than by errors in the observed frequencies. 
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